Coronal fuzziness modelled with pulse-heated multistranded loop systems 
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ABSTRACT 

Coronal active regions are observed to get fuzzier and fuzzier (i.e. more and more 
confused and uniform) in harder and harder energy bands or lines. We explain this 
evidence as due to the fine multi-temperature structure of coronal loops. To this end, 
we model bundles of loops made of thin strands, each heated by short and intense heat 
pulses. For simplicity, we assume that the heat pulses are all equal and triggered only 
once in each strand at a random time. The pulse intensity and cadence are selected 
so as to have steady active region loops (~ 3 MK), on the average. We compute the 
evolution of the confined heated plasma with a hydrodynamic loop model. We then 
compute the emission along each strand in several spectral lines, from cool (< 1 MK), 
to warm (2-3 MK) lines, detectable with Hinode/EIS, to hot X-ray lines. The strands 
are then put side-by-side to construct an active region loop bundle. We find that in the 
warm lines (2-3 MK) the loop emission fills all the available image surface. Therefore 
the emission appears quite uniform and it is difficult to resolve the single loops, while 
in the cool lines the loops are considerably more contrasted and the region is less fuzzy. 
The main reasons for this effect are that, during their evolution, i.e. pulse heating and 
slow cooling, each strand spends a relatively long time at temperatures around 2-3 
MK, and that it has a high emission measure during that phase, so the whole region 
appears more uniform or smudged. We make the prediction that the fuzziness should 
be reduced in the hot UV and X-ray lines. 

Subject headings: Sun: corona - Sun: X-rays 
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1. Introduction 

Coronal loops form the basic building blocks of the solar corona. The average temperature 
of these structures ranges from less than 1 MK in cool loops to more than 10 MK during flares 



(Vaianaetal. 



1973f) . The plasma inside coronal loops emits radiation mainly in the EUV and 
X energy band, with spectra depending on temperature. Hence, multiband and/or multiline 
observations are fundamental for the investigation of the thermal structure of the corona. 

Comparison of observations in different bands and different lines has shown that the confined 



corona, and in particular a ctive regions, appear to be fuzzier in hotter bands and lines, better 



defined in co oler ones (e.g. iBrickhou se & Sc hmel 



according to 



Tripathi et al. 



2006; 



Tripathi et al. 



2009). More specifically, 



(|2009|) . a loop bundle is fuzzier when it is difficult to resolve the 
single loops and less fuzzy when the separations between loops are more clear and with a higher 
contrast. This effect has been well observed with imaging instruments. 

o 

Images taken with the 284 A p assband filter of the Transition Region and Coronal 



Explorer (TRACE; 



Handy et al. 



Heliospheric Observatory (SOHO; 



filters sensitive to cooler plasma. 



1999), and of the 



Domingo et al 



UV Image Telescope (EIT) on the Solar and 



1995) appear to be fuzzier than in the other 



Brickhouse & Schmelzl (|2006 ) found that resonance scattering 



can be ruled out as the cause of this fuzziness. They also suggest that it may be due to the effect 



of filamentation, and particularly to the high emission measure 
preliminary modeling was performed for TRACE observation by 



between 2 and 3 MK. Some 



Patsourakos et al 



(2002) 



The Extreme-ultraviolet 



spacecraft (Kos ugi et al. 



framework (jTripathi et al 



maging Spectrometer (EIS; Culhanell2007[) . on board of Hinode 



2007), puts this evidence under a more quantitative and objective 



20091) . discovering that the same active region structures are clearly 



discernible in cooler lines ( ~ 1 MK), and are fuzzy at a higher temperature (~ 2 MK). 



The question we investigate here is whether the increasing fuzziness of coronal images in 
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hotter and hotter bands is really not an instrumental effect, but it is an intrinsic feature of the 
structure of coronal loop systems, and in particular to their fine transversal spatial and thermal 
structuring. In this sense, we will explore the hypothesis that the same loop system looks fuzzier 
and fuzzier at higher and higher temperatures, and not that this effect is due to two different classes 



Sakamoto et al 



of loo ps, one cooler and more clearly discernible, and the other hotter and fuzzier (I 
20091) . We will model coronal loops as composed by bundles of unresolved strands each ignited 
at random time by a rapid heat pulse, and compare their appearance in different spectral lines. 
Impulsively heat ed strands are commonly invoked in th e framework of the nanoflare coronal 



heating model (IParker 



1988 



Careill & Klimc huk 



2004J). The evolution of a single loop strand 



will be computed in detail by means of a numerical hydrodynamic loop model, and we will then 
form a loop system by collecting ensembles of strands. Each strand, at a given time, will be at a 
different phase of its evolution and will therefore emit differently. Globally the loop system will 
look different when observed in different spectral lines and we will study if the difference explains 
the different degree of fuzziness. 



In Sec. [2]we describe the model, in Sec. [3]the results are shown and discussed in Sec. HI 



2. The model 

We address a typical loop system of an active region. The system is a bundle of similar 
loops, each consisting of many thin strands. Each strand is pulse-heated. Before the heating, the 
strand is tenuous and cold, actually invisible. The loops consist of multitudes of strands, and, in 
principle, we should compute the hydrodynamic evolution for each of them. Instead of doing 
this, we make the simplifying assumption that the loops are made of similar strands and that 
each strand is heated by the same identical heat pulse, and therefore undergoes the same identical 
evolution. The heating rate averaged over the whole loop system is that of a steady-state loop at 3 
MK, i.e. a typical active region temperature. We assume a constant rate of heat pulses across the 
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loops. Therefore, the steady state is reached after a transient in which the number of strands with 
an on-going heat pulse increases linearly. The loop half-length is L = 3 x 10 9 cm. We assume that 
the loops stand vertical on the solar surface. 



The plasma confined in each strand transports energy and moves only along the magneti c 



1980; 


Peres et a 


1982 


, Doschek et al. 


1982; 


Nagai & Emslie 


1984; 


Fisher et a 


1985; 


Mac 


Meice 


1986; 


Gan et al. 


L991 


;] 


lansteenl 


1993; 


Betta et al. 


1997; 


Antiochos et al. 


1999 


;I 


Vluller et al. 


2003 



Bradshaw & Mason 



Betta et al 



1997L 



2003 



Bradshaw & Cargill 



2006), through the equations (|Peres et al 



19821 : 



dn dv 
dt ds 



(1) 



dv dp d ( dv 

nm H-J7 = —5- + nm H g+ — \p— 
dt os os \ os 



(2) 



ds dv „ (dv\ d I ^ sn dT 



(3) 



p =(1 +j3)nK B T 



(4) 



s = -p + nfix 



(5) 



w = -p + nfa 



(6) 



where n is the hydrogen number density; t is the time, s is the field line coordinate; v is 
the plasma velocity; m H is the mass of hydrogen atom; p is the pressure; g is the component of 
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gravity along the field line; n is the effective coefficient of compressional viscosity (including 
numerical viscosity); = njn is the ionization fraction where n e is the electron density; T is the 
temperature; k is the thermal conductivity 9 • 10~ 7 erg cm' 1 s' 1 K~ 7/2 ); K B is the Boltzmann 
constant; y is the hydrogen ioniza tion potential; P(T) are the radiative losses per unit emission 



measure ([Raymond & Smith 



1977h : Q (s, t) is the volumetric power input to the solar atmosphere: 



Q(s,t) =H steady + H f(t)g(s) (7) 

H steady is the steady heating term which balances radiative and conductive losses for the static 
initial atmosphere; the second term describes the heat pulse as a separable function of space and 
time. 

As mentioned above, the initial atmosphere of all the strands is tenuous and cool, since we 
want them to be virtually invisible in any relevant spectral band. On the other hand, the initial 
atmosphere has anyhow to sustain a high input energy; in particular, the chromosphere has to 
provide a mass amount large enough to support the evaporation driven by the heat pulse. This 
requires that the strand pressure cannot be too low. Apart from this, the fine details of the initial 
strand are not important because the energy input from the heat pulse overwhelms completely 
the initial energy budget and the evolution is largely independent of the initial atmosphere. A 
good compromise solution for this technical issue has been a loop strand with a base pressure of 
0.55 • 10 1 dyn cm -2 , which results in an apex temperature of 8.0 x 10 5 K. Since the emission of 
such a strand would be low but not be completely negligible in UV lines such as Mg VII and Si 
VI, for the sake of clarity we have decided anyway only for the analysis to set the emission of the 
in itial strands to zero at any wavelength. For the chromospheric part of the strand we use model F 



in 



Vernazza et al. 



(|1981 ). appropriate for active regions. The energy balance is strictly maintained 



at all times also in the chromosphere. 

The heat pulse is distributed uniformly along the loop strand: 
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1 





if < s < L 
elsewhere 



(8) 



The results do not change significantly if the he at pulses are deposited a t the loop footpoints, 



instead of uniformly, one of the options explored by 



The time dependence is a pulse function: 



Reale & Orlando (2008). 



fit) 



(9) 



1 if t < t < t + 60 sec 
any other t 

where t is the start time of heat pulse. The amplitude of the pulse is H = 0.38 erg cm~ 3 s' 1 . 
Note that the duration of heat pulse is short enough to have a multi-temperature loop system, and 
not short eno ugh to have large effects on the overall evolution and emission due to non-equilibrium 



of ionization (Reale & Orlando 



2008) 



The equation s are solved num erically by means of the Palermo-Harvard loop code 



(IPeres et al. 



both flaring (|Peres et al. 



982; 



Betta et al 



19871 : 



997) . This is a well tested and highly stable code, used for 



Betta et al 



Palermo-Harvard code has an adaptive grid (IBetta et al. 



2001) and quiescen t loops (IReale et al 



2000). The 



19971) . to better describe the steep 



gradients along the strands and during the evolution. For the ease of post-processing, the code 
output results were interpolated on a fixed equispaced grid for the post-processing. The grid is 
made of 1024 cells along the strand. 

The code output consists of temperature, density and velocity distributions along the loop 
strand sampled with a regular cadence during the strand evolution driven by the heat pulse. 

As mentioned above we model only one strand and replicate it so that one strand is different 
from the other only for the start time of the heat pulse, i.e. t . This choice minimizes the number 
of free parameters, and a single simulation was enough for our further analysis. The simulation 
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computes the evolution of the strand for 2000 s after the start of the heat pulse, and the solutions 
are sampled every 1 s. 

We have also assumed that all the strands are strictly out of phase, i.e. only one strand is 
switched on at any time. At time t = all strands are switched off. Then the heat pulses gradually 
turn them on, one after the other. We have chosen to turn them on with constant cadence to have a 
smooth transition from no bright strands to maximum number of bright strands (Fig. [4]). 

At a certain time, part of the strands will still be off, in a few the heat pulse will be on, the 
others will be cooling after the heat pulse has ended. Since all strands have the same evolution, 
at each time each strand will be at a different phase of the same evolution. A whole loop at a 
given time is therefore simulated as a set of otherwise identically evolving strands, but each with 
a different heating start, all glued together. This approach is not completely new; outputs of the 
same loop simula t ion have been put togethe r in the past to simulate m ultistrand nanoflaring loops 



(IPeres et al 



1993 



Warren et al 



2002, 



2003 



Winebarger et al 



2003allb|) but in all cases they were 



averaged to obtain the effective aspect and evolution of a loop as a whole. The novelty of our 
approach is that we do not model each loop as whole, but we keep the information of its spatial 
longitudinal and transversal distribution. 

In our model one heat pulse is switched on every one second in one new strand. One loop 
consists of « 60 strands and we put « 30 loops one by the other. Every second there are 60 
strands where the heat pulse is on out of a total of 2000. In the steady-state the loop heating per 
unit volume is the refore H = H Q x 60/2000 ~ 0.011 erg cm 3 s _1 . According to scaling law 



(Rosner et al. 



1978[), this steady heating sustains a loop with an apex temperature of « 3.8 MK 



corresponding to an average temperature along the loop of about 3 MK(Fig. HI). 

When we put all strands side by side, we obtain a mesh 1024 x 2000 grid cells. The 
transversal width of each strand is constant along the loops. Therefore it is only a multiplicative 
constant which determines the cross-section of the entire loop system. 
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From model results we synthesized the loop emission in different spectral lines. We selected 
the spectral lines, in Tab. \T\ with a peak emis sivity in a wide rang e of temperatures. A subset 



(the first seven) are detected with Hinode/EIS (ITripathi et al. 



(|2009|) ). which provide our reference 



observational evidence. We also include hotter lines detectable in the UV (e.g. Fe XFX, Fe XX, 
Fe XXIII from AIA on board the SDO mission) and X-ray band. 



The emissivity in these lines (G(T), Fig.Q]) is take n from the CHIA NTI database ( Dere et al 



1997 



Young et al 



2003|), with coronal abundances of 



10 9 cm 3 ), and ionization fraction of iMazzotta et al. 
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Feldmanl (|1992[) (assuming a density of 



FeXV 



CO 

E 
o 

D) 

! 

CD 



CO 



o 

CO 
CO 



10° 
10" 2 
10" 4 
10" 6 
10" 8 



Fe-XVJ/ 



E 10 - 10 

LJJ 1 w 



10 



12 




FeXXV 



FeXX'' 



-feXXUl 
■CaXIX 

':;;-.'.';;--.,sxv 

"'•■•'•'.':;SiXII 

MgX 



FeXIX 



5.0 5.5 6.0 



6.5 7.0 
Log T[K] 



7.5 8.0 8.5 



Fig. 1. — Profile of emission per unit of emission measure vs log T\K] for the va rious spectral 



lines. We make a distinction between lines for comparison with ITripathi et al 



(2009) (solid lines), 



lines in AIA filter bands on board of SDO mission (dashed lines), and hot X ray lines (dotted lines) 



The emission from each grid cell is computed as: 
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I(T,n e ) = n]G{T)dV (10) 

*J V pixel 

where V is the volume and V pixe i is the volume enclosed by a single pixel. 

To obtain a final image to be compared with observational data we follow the conceptual 
steps shown in Fig.[2l We first compute the emission along each strand in a given line. The strand 
emissions can be put in the form of a 1-D image consisting of a strip of pixels on an appropriate 
grey scale. We obtain 2000 pixel strips. We then put them side by side and obtain a 1024 x 2000 
pixel image. We have chosen to rebin our images by summing over bundles of 65 strands, so to 
obtain a collection of about 30 parallel loops. 



3. The results 



T he evolution of nanoflaring plasma confined in coronal loop is well known from previou s 



works ( Peres et al. 



1993 



Warren et al. 



2002, 



2003; 



Patsourakos & Klimchuk 



2005 



200 51) and is simi 



Peres et al. 



dl982h : 



a r to that of proper flari ng loops, although on a smaller scale (e.g., 



Testa et al. 



Nagai ( 



1980); 



Reale & Peresl (|1995|) ). We show the evolution of the density and temperature 
along a single strand in Fig. [3l The temperature soon settles to about 10 MK along most of the 
strand, due to the strong heat pulse; then it slowly decreases as expected with an e-folding time 
scale given by t s « 500 Lg/ y[Tl « 500 s (Reale 2007, 2009) where T 6 and L 9 are the maximum 
temperature and the length of the strand measured in units of 10 6 K (MK) and 10 9 cm respectively. 
At t = 2000 s the strand has cooled below the temperature it had before the heat pulse. At t = 20 s 
the density plot clearly shows a strong evaporation front coming up from the chromosphere (the 
density of the front jumps by more than a factor 10). At t = 100 s the front has filled the whole 
strand with some extra accumulation in the apex region. At t = 400 s the density distribution has a 



shape similar to the initial equilibrium one, but settled around a maximum value of ~ 10 cm . 
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J, 





Fig. 2. — Flowchart showing how an loop bundle image is obtained from the emission of each individual 
loop in a given spectral line. We first compute the spatial distribution of the emission along each strand, we 
then put the emission in the form of an image made of a pixel strip, we put all pixel strips side by side; we 
then reconstruct the aspect of the whole loop system, then we group the strands into loops. Finally we bend 
the image to obtain a loop-like shape. 
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Then the density begins to decrease and at the final time computed it is lower, by a factor ~ 5 than 
the maximum, still much higher than the initial value. 

As described in Sec.[2l this evolution is replicated in all strands with different relative timing, 
depending on the time the heat pulse is switched on. As mentioned above, for a realistic situation, 
we have simulated the ignition of the whole loop system, with a gradually increasing number of 
simultaneously heated strands. The transient evolution will be the subject of a future work, but 
Fig. H] shows the initial evolution of the temperature averaged over the whole loop system. In the 
following we will focus our attention to the final time (t = 2000 s) when the loop system enters a 
presumably long steady state, in which the heat pulse 'storm' (i.e. the 2000 pulses in total) repeats 
continuously. In a way, we are therefore implicitly assuming that the time taken to re-energise the 
magnetic field (e.g. by twisting, braiding) is about 2000 s. 

Fig. [5] shows the distribution of the emission mea sure of the en tire collection of strands, 



Cargilll 



1994|) at this time. We show the 



i.e. of the whole loop system, versus temperature (e.g. 
distribution of the coronal part only, i.e. upper 90% of the loop bundle, excluding the lower layers. 
The peak of the distribution is around 3 MK, as planned. We also see that the distribution is quite 
broad; the flatter tail is toward the cool part, but there are also significant hot components up to 
10 MK, as expected, due to the presence of the strong heat pulses. These hot components are of 
course minor, due to the small duty cycle of the heating phase with respect to the whole evolution 
of the strand. 

Using the method described in the Sec [2l we synthesize maps of emission for all the lines 
considered. Fig. [6](left column) shows a subset of them, namely the map in cool EIS lines, i.e. Mg 
VH A278 (log T[K] = 5.8), medium temperature EIS line Fe X A1S6 (log T[K] = 6.0), warm 
EIS line Fe XV A2S4 (log T[K] = 6.4), and hot X-ray line Ca XIX ,13.21 (log T[K] = 7.4). 



The images on the 



Tripathi et al 



eft column are analogous, and can be compared, to observed ones (e.g. 



(|2009|) ). The grey scale is logarithmic and spans a factor 10 in all maps and plots. 
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Fig. 3. — Temperature (a), and density (b) along half of a single strand at times t = 0, 20, 100, 400, 
800, 2000 s (thinner and thinner lines with progressing time). The heat pulse starts at time t = 
and lasts 60 s. 
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500 1000 1500 2000 

Time (s) 

Fig. 4. — Evolution of the loop average temperature from off-state to steady-state. 




Fig. 5. — Distribution of emission measure vs Temperature at t = 2000 s. We show the distribution 
of the coronal part only, i.e. the upper 90% of the loops. 
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(g) Fe XXIII (h) Fe XXIII 

Fig. 6. — Loop bundle emission in different (labeled) spectral lines. In the images (left column) we 
show the upper 90% of the loop. The grey scale is logarithmic, and covers a factor 10 of intensity; 
black corresponds to the maximum of emission and white to the minimum. The plots in the right 
column are the emission profiles along the lines marked in the images of the left column, in the 
same logarithmic scale. 
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This is a customary choice for showing observed data (e.g. iTripathi et al.1 (|2009|)). In the first three 
lines, the emission decreases from the base to the top of the loop, due to the density decreasing as 
well. We show the intensity of the upper 90% of the loop system, i.e. the coronal part only. In this 
way we exclude the moss emission from cromosphere and transition region, which, nevertheless, 
would appear saturated in this color scale. This may indicate that the related emission measure 
predicted by the model is too high, which i s a w ell known effect of 1-D hydrodynamic loop 



modeling (e.g. 



Warren & Winebarger 



2006 



2007|), an d ma y point to the need for a more accurate 



description of the low loop atmosphere. This is however not required in this context. In the hot Fe 
XXIII line the opposite occurs, i.e. the brightness increases from the base to the top of the loop 
system, because the line is sensitive to plasma with log T > 6.9, which is practically never found 
in the low part of the loop system. For our specific question, from the comparison of the first three 
lines we can immediately see that we can easily resolve the single loops in the cool line image, 
that this is more difficult in the second image, and that the Fe XV image is much more uniform 
and does not practically allow us to resolve the structures. Interestingly, the fourth image in the 
hot UV line shows again very well contrasted structures. This occurs also for the other hot lines, 
although with different levels of contrast. The plots in the right columns of Fig. |6]quantify better 
this effect: the profile along the cut in the cool line shows large fluctuations, with large amplitude 
and high spatial frequency, in the warm line pro file the fluctuations a re of smaller amplitude. 



According to the operative definition reported in 



Tripathi et al 



(|2009r) . the loop system in the Fe 



XV image is fuzzy, the loop syst em in the Mg VII im age is not. This is in total agreement with the 



observational evidence shown in 



Tripathi et al. 



(120091) 



We have also devised a figure for fuzziness to be computed in the cross-section line: 



1 



(11) 



where /, is the emission in the spectral line in the i-th pixel along the cross-section, / is the 
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Fig. 7. — Scatter plot of the values of density and temperature, obtained from averaging on bins 
of 32 grid points in the whole model strand evolution (upper 90% part of the loops). The dashed 
line is the locus of density and temperature at t he apex of hydrosta tic loops, with half-length 



L = 3.0 x 10 9 cm, according to the scaling laws of iRosner et al. 



(1978). 
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average value for and n is the number of pixel in the region marked by the line (left panels in 
Fig. [6]). Clearly, the larger this quantity, the more contrasted and well resolved are the loops, and 
the smaller is the fuzziness. Table \T\ includes the values of this figure for all the spectral lines 
analyzed. All these value are calculated along the same cut shown in Fig. [6l 

The values are quite high for the lines at low temperatures (logT < 6; Mg VII, Si VII, Fe 
X), very high for the lines at very high temperature (logT > 6.9; e.g. S XV, Ca XIX, Fe XXIII, 
Fe XXV), lower for the warm lines (around 3 MK; especially FeXV, FeXVI, Ne IX). We have 
applied the same analysis for different binnings, and still find similar results. 

The reason why the loop system appears more uniform at warm temperatures is clear from 
inspection of Fig. [7] In the figure each data point marks a value of density and temperature 
obtained from averaging over sections of 32 grid points along a strand (upper 90%) in each output. 
The figure clearly shows that around 3 MK the density of data points is higher, and at higher 
values of density. The implication is that each strand spends a long time, and with a high emission 
measure, at an average temperature of about 3 MK, and therefore at this temperature more strands 
appear on the average brighter, i.e. the loop system is more uniformly bright. The plasma does 
not stay long at a higher temperature, i.e. around the heat pulse, and moreover during this phase, 
the emission measure is not high because the evaporation is still in an early phase. The plasma 
does instead stay for a long time at lower temperature, but the emission measure becomes very 
low, not contributing much to the overall emission. Fig. [7] also shows that this model fits well the 
obser vational constra int of alternating cool/overdense - hot/underdense status of coronal loops 



(e.g. 



Klimchuk 



2006). 



It is easy to imagine extrapolating this result to whole active regions, threaded with many 
thousands upon thousands of magnetic field lines. There are fewer field lines populated with 
low and high temperature plasma and so distinct loop structures appear. There are many more 
field lines populated with plasma in the 2 - 3 MK range and so fuzzy emission appearing to be 
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Table 1 : Variation of emission in the lines analyzed. 



Line 


Wavelength (A) 


Emission temperature (logT[K]) 


(Till 


Mg VII 


278 


5.8 


0.28 


Si VII 


275 


5.8 


0.26 


FeX 


186 


6.0 


0.23 


FeXII 


195 


6.1 


0.22 


Fe XIII 


196 


6.2 


0.21 


Fe XIV 


264 


6.3 


0.20 


Fe XV 


284 


6.4 


0.18 


Fe XVI 


361 


6.5 


0.15 


NelX 


13.4 


6.6 


0.12 


Mg At 


n o 1 
9.31 


6.8 


A 1 A 
(J. 14 


Fe XIX 


133 


6.9 


0.26 


FeXX 


133 


7.0 


0.39 


Si XIII 


6.69 


7.0 


0.21 


Fe XXIII 


133 


7.1 


0.55 


S XV 


5.10 


7.2 


0.30 


Ca XIX 


3.21 


7.4 


0.46 


FeXXV 


1.87 


7.6 


1.0 
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composed of unresolved structures would arise as observed. 



4. Discussion and Conclusions 

In this work we investigate the reason why coronal loop systems, or entire active regions, 
look fuzzier in warm ~ 2 - 3 MK spectral lines than in cooler < 1 MK lines. This evidence has 
been well known since the first X-ray and UV missions, e.g. Skylab, but it has been rece ntly put 
on a more established ground from Hinode/EIS spectral observation (|Tripathi et al.ll2009|) . Our 
basic scenario is that coronal loops consist of bundles of thin strands, each of thickness below 
the instrumental spatial resolution, and that each strand is heated up to about 10 MK by a strong 
and fast heat pulse, i.e. the loops are heated by a storm of nanoflares. The plasma is confined 
in each strand, so that it evolves as an independent atmos phere, and can be modeled with loop 
hydrodynamics (see also IPatsourakos & Klimchuk! (|2007l) for a conceptually similar approach). 
Our choice has been to assume that the strands are all heated once and by the same heat pulse, 
lasting 60 s, occurring at a different random time for each strand, with a cadence and an intensity 
adequate to maintain the loop at ~ 3 MK on average. This model fits well the observational 
constraint of hot/underdense-cool/overdense cycles. We have then collected 2000 different strands 
to form a loop system, and derived synthetic images of the loop system when it reaches steady 
state in several relevant spectral lines. In our opinion, the images synthesized from our model 
unequivocally show the same "fuzziness" in the same warm lines as observed with EIS, and the 
same better definition in the cool lines as observed with EIS. In other words, our model is able 
to explain the evidence. Of course, it explains also the effect as observed in narrow-band XUV 
instruments such as the normal-incidence imaging telescopes, TRACE and SoHO/EIT. We have 
also provided quantitative figures to this effect. The basic reason why this model works is that, 
in spite of the short heat pulse, the strands spend a long time with a high emission measure at a 
temperature around 3 MK, much less time when plasma is hotter and long time, but with much 



-22- 



less emission measure when the plasma is cooler. 



So the loop systems appear more uniform aro und 3 MK, and this higher filling factor gives 



the impression of "fuzziness", as described in 



Tripathi et al 



(2009). In cooler lines we are able to 



resolve better the loops, which appear more contrasted and with better defined boundaries. 



This "fuzziness" is different from the one intended in lSakamoto et al.l (|2009i) : we address 
the question why the same loop s appear diff erent in different line s, which is exactly the evidence 



reported in 



Tripathi et al 



(2009), whereas 



Sakamoto et al 



(2009) address the evidence that hot 



loops appear fuzzier than cooler, not co-spatial, and therefore different, loops. 

There are several limitations in our approach. Small duration heat pulses are a necessary 
ingredient to have a multi-temperature loop system. Here we assume that the heat pulses are the 
same in all strands, only their timing is different, and that, at regime, they occur with a constant 
time average. We expect some broadness in the distribution of the heat pulses, in duration, 
intensity and average cadence. Other simplifications concern other assumptions on the heat 
pulses. We are assuming a very simple time-shape of the heat pulse, with abrupt on and off. 
This should not affect the results, because we analyze the strands evolution on a much longer 
time scale. The heat puls e is long enough not to have significative effects due to delay to reach 



equilibrium of ionization (IReale & Orlando 



2008). We also assume that each strand is heated only 



once, i.e. there is no reheating of the same strand. This is a realistic assumption in the framework 
of progressive magnetic reconnection occurring in gradually twisting loops, due to footpoint 
photospheric motions. Presumably there can only be one heating event per strand since it is 
destroyed after magnetic reconnection. Nevertheless, note that the reconnecting strands form two 
new (relaxed) strands which are then presumably subject to the same twisting/braiding motions 
which originally energised them, leading to eventual reconnection and another heating event. In 
principle, this process could be repeated many times for a single pair of strands. The magnetic 
flux in this case is not destroyed, simply reconfigured. 



The hydrodynamic description of the confined plasma should hold quite firmly in the 

ow loop section 



Gabriel 



1976), but 



individual loop strands, as in monolithic loops. Maybe the description of the 
should, for instance, include the possible tapering in the transition region (e.g. 
differences are expected only in the very low part of the loops and, therefore, should not affect the 
results presented here. 

Our analysis provides also an interesting prediction: we expect high contrast to show in very 
hot lines, such as the Fe XXIII line and the Ca XIX line (typical flare lines), even stronger than 
in the cool lines. Therefore, forthcoming SDO observations and desirable future X-ray spectral 
observations of the quiet Sun, and in particular of active regions, at high sensitivity and enough 
spatial resolution, may provide an important confirmation of our scenario. 

In this analysis, we have studied a coronal loop at steady-state, i.e. when its light curve and 
temperature is constant on average. However, our model is useful also to study the time evolution 
of the multistranded loops, and we will further develop our analysis in this perspective as the next 
step. 
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